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The production of quarkonium in heavy ion collisions is studied at RHIC and LHC 
energies. General arguments indicate that, due to initial production of multiple quark 
pairs in each central collision, the final population of quarkonium may exhibit significant 
enhancements over straightforward extrapolation of behavior at SPS energy. Explicit 
| calculations based on both a statistical hadronization picture and a kinetic formation 
mechanism in a deconfined state verify these general expectations. Such enhancements will 
alter the nature of how quarkonium yields may be used as a signature of deconfinement. 



1. GENERIC EXPECTATIONS FOR LARGE N r , 



At heavy ion collision energies at RHIC and LHC, the initial number of heavy quark 
pairs produced in each collision will be qualitatively different than the number produced 
at SPS energies. For example, the number of charm quark pairs is expected to be of 
^ the order of ten at RHIC and several hundred at LHC in the most central collisions [|TJ. 
Let us attempt to extract features of J/ip formation from the initially-produced charm 
quarks which are independent of detailed dynamics. We consider scenarios in which the 
formation of J/ip is allowed to proceed through any combination of one of the N c charm 
quarks with one of the N 5 anticharm quarks which result from the initial production of 
N c5 pairs in a central heavy ion collision. For a given charm quark, one expects then that 
the probability V to form a J/ip is just proportional to the number of available anticharm 
quarks relative to the number of light antiquarks, 

V oc N- c /N^- s « N cd /N ch , (1) 

where we normalize the number of light antiquarks by the number of produced charged 
hadrons. Since this ratio is generally very small, one can simply multiply by the number 
of available charm quarks N c to obtain the total number of J/ip expected in a given event. 

N JH , cx Nj/N ch} (2) 

where the use of the initial values N c5 = N c = N 5 is again justified by the relatively small 
number of bound states formed. For an ensemble of events, the average number of J / ip 
per event is calculated from the average value of initial charm < N c5 >, and we neglect 
fluctuations in N^. 

< J/-0 >= A(< N c5 > +1) < N cS > /N ch , (3) 
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where we place all dynamical dependence in the parameter A. [One can extend this 
formula to the case where J /if) formation is effective not over the entire rapidity range 
Ytotah but only if the quark and antiquark are within the same rapidity interval Ay. In 
this case one makes the replacement < N c£ > +1 — >< N cS > +Y tota i / Ay .] 

The essential property of this result is that the growth with energy of the quadratic 
dependence on total charm M is expected to be much stronger than the corresponding 
growth of total particle production in heavy ion collisions @. J /if) production without 
this quadratic mechanism is typically some small energy-independent fraction of total 
initial charm production 0, so that we can expect the quadratic formation to become 
dominant at high energy. 

We show numerical results in Table 1 for these quantities with a prefactor A of unity. 
Estimates for the charm and particle numbers are very approximate, but serve to show 
the anticipated trend with energy. At SPS, this formation mechanism is most probably 
insignificant. At RHIC it is comparable with "normal" formation, while at LHC one 
might expect it to be dominant. Of course, the exact result will depend on the details of 
the physics which controls the formation. 



Table 1 



Comparison of J /if) formation 


variation with energy 






SPS 


RHIC 


LHC 


(GeV) 


18 


200 


5500 


< N c5 > 


0.2 


10 


200 




1350 


3250 


16500 


< N JH> > 


0.00018 


0.034 


2.4 


^initial 


0.0012 


0.06 


1.2 



2. STATISTICAL HADRONIZATION MODEL 

This model is motivated by the success of attempts to explain the relative abundances 
of light hadrons produced in high energy interactions in terms of the predictions of a 
hadron gas in chemical and thermal equilibrium Such fits, however, are not able to 
describe the abundances of hadrons containing charm quarks. This can be understood in 
terms of the long time scales required to approach chemical equilibrium for heavy quarks. 
However, it is expected that for high energy heavy ion collisions the initial production of 
charm quark pairs exceeds the number expected at chemical equilibrium as determined 
by the light hadron abundances. The statistical hadronization model || assumes that 
at hadronization these charm quarks are distributed into hadrons according to chemical 
equilibrium, but adjusted by a factor 7 C which accounts for oversaturation of charm. 
One power of this factor multiplies a given thermal hadron population for each charm 
or anticharm quark contained in the hadron. Thus the relative abundance of J /if) to 
that of D mesons, for example, will be enhanced in this model. The enhancement factor 
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is determined by conservation of charm, again using the time scale argument to justify 
neglecting pair production or annihilation before hadronization. 

N c5 = \ lc N * l{ -? cN J en \ + lc 2 N hldden , (4) 

where N open and N hidden are calculated in the thermal equilibrium model grand canonical 
ensemble. The ratio of Bessel functions is an approximate correction factor |J which 
converts to canonical ensemble particle numbers. This factor will only be significant for 
non-central collisions or lower energies in which the absolute number of charm pairs is 
quite small f7|. 

We show in Figure [I] calculated 7 C values. At large total charm there is linear behavior, 
which leads to the quadratic growth of J/ip formation. It is interesting to note that the 
canonical suppression effect can be very well approximated by using a grand canonical 
formalism, supplemented by the averaging process for N% 5 over a distribution of events. 

Figure |2| shows several applications for RHIC conditions. The centrality dependence 
is modeled by the number of nucleon participants, and one sees the change in shape 
due to the transition between canonical and grand canonical formalism. The absolute 
magnitudes of < J/ip >/charm are comparable with the initial production estimates of a 
fraction of a percent, indicating that this process may overwhelm suppression for central 
collisions at RHIC. 




Figure 1. Relation between charm enhance- Figure 2. Ratio < J/ip > over initial charm 
ment factor and total number of charm quarks at RHIC for several applications of the sta- 
for statistical hadronization model. tistical hadronization model. 



3. KINETIC FORMATION MODEL 



In this model, we investigate the possibility to form J/ip directly in a deconfined 
medium. The formation will take advantage of the mobility of initially-produced charm 
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quarks in a spatial region of deconfinement, such that all combinations of a charm plus 
anticharm are allowed to form a bound state. For the purposes of this study, we consider 
a physical picture of deconfinement in which quarkonium is suppressed via collisions with 
free gluons |J. The dominant formation process, in which a quark and an antiquark in 
a relative color octet state are captured into a color singlet bound quarkonium state and 
emit a color octet gluon, is simply the inverse of the breakup reaction. It is then an 
inevitable consequence of this picture of suppression that the corresponding formation 
process must also take place. 

The time evolution if the J /if) population is given by the rate equation 

= \ F N c N s [V(r)]- 1 - X D N J/fPg , (5) 

with p g the number density of gluons. The reactivity A is the reaction rate (av Te i) averaged 
over the momentum distribution of the initial participants, i.e. c and c for Xp and J/if> 
and g for Ad- The gluon density is determined by the equilibrium value in the QGP at 
each temperature. The system undergoes a longitudinal isentropic expansion, which fixes 
the time-dependence of the volume V(r). It is evident that the solution of Equation [5] 
grows quadratically with initial charm N c5 , as long as the total J/if> <C N c5 . 

Nj/Arf) = < r f) x Wj/^ro) + N 2 C - C P A F [V(r) e(r)]- 1 dr], (6) 

J TO 

where Tf is the hadronization time determined by the initial temperature (T is a variable 
parameter) and final temperature (Tf = 150 MeV ends the deconfining phase). The 

function e(ry) = e ^ T o x ° p3dT would be the suppression factor in this scenario if the 
formation mechanism were neglected. 

The momentum distribution of the charm quarks is allowed to vary over a wide range 
of possibilities. At one extreme we use a thermal equilibrium distribution at the QGP 
temperature. We also consider a distribution unchanged from that introduced in the 
initial perturbative QCD processes, plus several intermediate distributions with decreasing 
rapidity widths. The centrality dependence can also be modeled, using nuclear collision 
geometry to set the various spatial scales. For details, see Ref. [|9|,|T0|1 . 

We show the resulting ratios Nj/^/N c5 in Figure |3| (for RHIC) and Figure £| (for LHC) 
for a full range of the initial charm quark rapidity distributions. Shown for comparison 
are some statistical hadronization model results, plus an indication of the value if only 
initial (vacuum) production were present. It is evident that although the absolute values 
depend strongly on the charm momentum distribution, the magnitudes are substantially 
above that expected if only initial production plus plasma suppression were present. The 
centrality dependence is quite striking. Even after normalizing to initial charm, the ratio 
rises with increasing centrality. 



4. SUMMARY 

Expectations based on general grounds for enhanced formation of heavy quarkonium in 
relativistic heavy ion collisions have been verified in two different models. In particular, 
one expects at RHIC and LHC to see an enhancement in the heavy quarkonium formation 
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Figure 3. Predictions of kinetic model for 
< J/ip > over initial charm at RHIC energy. 



Figure 4. Predictions of kinetic model for 
< J/ip > over initial charm at LHC energy. 



rate, even when compared to unsuppressed production via elementary nucleon-nucleon 
collisions in vacuum. The magnitude of this effect is expected to grow with the centrality 
of the heavy ion collision, just opposite to the predictions of various suppression scenarios. 
The physics bases for these models, however, are quite distinct. Their differences should 
manifest themselves in details of the magnitudes and centrality dependence. In this 
regard, it is essential to have a simultaneous measurement of open flavor production to 
serve as an unambiguous baseline. 

REFERENCES 

1. P. L. McGaughey, E. Quack, P. V. Ruuskanen, R. Vogt, and X.-N. Wang, in "Hard 
Processes in Hadronic Interactions", Int. J. Mod. Phys. A10 (1995) 2999. 

2. X.-N. Wang and M. Gyulassy, Phys. Rev. Lett. 86 (2001) 3496. 

3. R. Gavai, D. Kharzeev, H. Satz, G. Schuler, K. Sridhar, and R. Vogt, in "Hard 
Processes in Hadronic Interactions", Int. J. Mod. Phys. A10 (1995)3043. 

4. See for example P. Braun-Munzinger, D. Magestro, K. Redlich, and J. Stachel, Phys. 
Lett. B518 (2001) 41, and references therein. 

5. P. Braun-Munzinger and J. Stachel, Phys. Lett. B490 (2000) 196. 

6. J. Cleymans, K. Redlich, and E. Suhonen Z. Phys. C51 (1991) 137. 

7. M.I. Gorenstein, A.P. Kostyuk, H. Stocker, and W. Greiner, Phys. Lett. B509 (2001) 
277; P. Braun-Munzinger and J. Stachel, Nucl. Phys. A690 (2001) 119. 

8. D. Kharzeev and H. Satz, Phys. Lett. B334 (1994) 155. 

9. R. L. Thews, M. Schroedter, and J. Rafelski, Phys. Rev. C63 (2001) 054905. 

10. R. L. Thews and J. Rafelski, J/ij) Production at RHIC in a QGP, |hep-ph/0104025| , 
to be published in proceedings of QM2001. 



